Epigenetic features, methods, and implementations associated with COVID-19

The infection and life cycle of severe acute respiratory syndrome coronavirus 2 are widely studied, yet multiple gaps exist in the knowledge that affects therapeutic developments against coronavirus disease 2019 (COVID-19). Predominantly caused by a respiratory virus, COVID-19 is not restricted to the respiratory tract but affects multiple organs of the body including the cardiovascular, neurological, immunological, and renal systems. COVID-19 affects all age groups, although the elderly population inherently presenting with multiple comorbidities are disproportionately affected. The majority of the patients experience mild symptoms, although moderate, severe, and critical symptoms occur in a smaller group of patients. Interestingly, the effects of the disease can be acute or chronic and present an ongoing health care challenge. Epigenetic mechanisms of COVID-19 (DNA methylation, histone posttranslational modifications, histone citrullination, etc.) are an emerging field and present enormous potential toward the medical management of COVID-19. Angiotensin converting enzyme 2, an important protein in the cardiovascular system, is a receptor for viral entry into cells, and the epigenetic processes that regulate this protein have been widely studied. Identification of the epitranscriptomic profile has led to the identification of putative biomarkers for disease diagnosis and trials of novel epidrugs for targeted therapy. © 2023 Elsevier Inc. All rights reserved.

No evidence for epitranscriptomic m5C modification of SARS-CoV-2, HIV and MLV viral RNA.

The addition of chemical groups to cellular RNA to modulate RNA fate and/or function is summarized under the term epitranscriptomic modification. More than 170 different modifications have been identified on cellular RNA, such as tRNA, rRNA and, to a lesser extent, on other RNA types. Recently, epitranscriptomic modification of viral RNA has received considerable attention as a possible additional mechanism regulating virus infection and replication. N6-methyladenosine (m6A) and C5-methylcytosine (m5C) have been most broadly studied in different RNA viruses. Various studies, however, reported varying results with regard to number and extent of the modification. Here we investigated the m5C methylome of SARS-CoV-2, and we reexamined reported m5C sites in HIV and MLV. Using a rigorous bisulfite-sequencing protocol and stringent data analysis, we found no evidence for the presence of m5C in these viruses. The data emphasize the necessity for optimizing experimental conditions and bioinformatic data analysis.

Epitranscriptome: Review of Top 25 Most-Studied RNA Modifications.

The alphabet of building blocks for RNA molecules is much larger than the standard four nucleotides. The diversity is achieved by the post-transcriptional biochemical modification of these nucleotides into distinct chemical entities that are structurally and functionally different from their unmodified counterparts. Some of these modifications are constituent and critical for RNA functions, while others serve as dynamic markings to regulate the fate of specific RNA molecules. Together, these modifications form the epitranscriptome, an essential layer of cellular biochemistry. As of the time of writing this review, more than 300 distinct RNA modifications from all three life domains have been identified. However, only a few of the most well-established modifications are included in most reviews on this topic. To provide a complete overview of the current state of research on the epitranscriptome, we analyzed the extent of the available information for all known RNA modifications. We selected 25 modifications to describe in detail. Summarizing our findings, we describe the current status of research on most RNA modifications and identify further developments in this field.

The epitranscriptome of Vero cells infected with SARS-CoV-2 assessed by direct RNA sequencing reveals m6A pattern changes and DRACH motif biases in viral and cellular RNAs.

The epitranscriptomics of the SARS-CoV-2 infected cell reveals its response to viral replication. Among various types of RNA nucleotide modifications, the m6A is the most common and is involved in several crucial processes of RNA intracellular location, maturation, half-life and translatability. This epitranscriptome contains a mixture of viral RNAs and cellular transcripts. In a previous study we presented the analysis of the SARS-CoV-2 RNA m6A methylation based on direct RNA sequencing and characterized DRACH motif mutations in different viral lineages. Here we present the analysis of the m6A transcript methylation of Vero cells (derived from African Green Monkeys) and Calu-3 cells (human) upon infection by SARS-CoV-2 using direct RNA sequencing data. Analysis of these data by nonparametric statistics and two computational methods (m6anet and EpiNano) show that m6A levels are higher in RNAs of infected cells. Functional enrichment analysis reveals increased m6A methylation of transcripts involved in translation, peptide and amine metabolism. This analysis allowed the identification of differentially methylated transcripts and m6A unique sites in the infected cell transcripts. Results here presented indicate that the cell response to viral infection not only changes the levels of mRNAs, as previously shown, but also its epitranscriptional pattern. Also, transcriptome-wide analysis shows strong nucleotide biases in DRACH motifs of cellular transcripts, both in Vero and Calu-3 cells, which use the signature GGACU whereas in viral RNAs the signature is GAACU. We hypothesize that the differences of DRACH motif biases, might force the convergent evolution of the viral genome resulting in better adaptation to target sequence preferences of writer, reader and eraser enzymes. To our knowledge, this is the first report on m6A epitranscriptome of the SARS-CoV-2 infected Vero cells by direct RNA sequencing, which is the sensu stricto RNA-seq.

Epitranscriptomics of SARS-CoV-2 Infection.

Recent studies on the epitranscriptomic code of SARS-CoV-2 infection have discovered various RNA modifications, such as N6-methyladenosine (m6A), pseudouridine (Ψ), and 2'-O-methylation (Nm). The effects of RNA methylation on SARS-CoV-2 replication and the enzymes involved in this mechanism are emerging. In this review, we summarize the advances in this emerging field and discuss the role of various players such as readers, writers, and erasers in m6A RNA methylation, the role of pseudouridine synthase one and seven in epitranscriptomic modification Ψ, an isomer of uridine, and role of nsp16/nsp10 heterodimer in 2'-O-methylation of the ribose sugar of the first nucleotide of SARS-CoV-2 mRNA. We also discuss RNA expression levels of various enzymes involved in RNA modifications in blood cells of SARS-CoV-2 infected individuals and their impact on host mRNA modification. In conclusion, these observations will facilitate the development of novel strategies and therapeutics for targeting RNA modification of SARS-CoV-2 RNA to control SARS-CoV-2 infection.

RESIC: A Tool for Comprehensive Adenosine to Inosine RNA Editing Site Identification and Classification.

Adenosine to inosine (A-to-I) RNA editing, the most prevalent type of RNA editing in metazoans, is carried out by adenosine deaminases (ADARs) in double-stranded RNA regions. Several computational approaches have been recently developed to identify A-to-I RNA editing sites from sequencing data, each addressing a particular issue. Here, we present RNA Editing Sites Identification and Classification (RESIC), an efficient pipeline that combines several approaches for the detection and classification of RNA editing sites. The pipeline can be used for all organisms and can use any number of RNA-sequencing datasets as input. RESIC provides (1) the detection of editing sites in both repetitive and non-repetitive genomic regions; (2) the identification of hyper-edited regions; and (3) optional exclusion of polymorphism sites to increase reliability, based on DNA, and ADAR-mutant RNA sequencing datasets, or SNP databases. We demonstrate the utility of RESIC by applying it to human, successfully overlapping and extending the list of known putative editing sites. We further tested changes in the patterns of A-to-I RNA editing, and RNA abundance of ADAR enzymes, following SARS-CoV-2 infection in human cell lines. Our results suggest that upon SARS-CoV-2 infection, compared to mock, the number of hyper editing sites is increased, and in agreement, the activity of ADAR1, which catalyzes hyper-editing, is enhanced. These results imply the involvement of A-to-I RNA editing in conceiving the unpredicted phenotype of COVID-19 disease. RESIC code is open-source and is easily extendable.

The Architecture of SARS-CoV-2 Transcriptome.

SARS-CoV-2 is a betacoronavirus responsible for the COVID-19 pandemic. Although the SARS-CoV-2 genome was reported recently, its transcriptomic architecture is unknown. Utilizing two complementary sequencing techniques, we present a high-resolution map of the SARS-CoV-2 transcriptome and epitranscriptome. DNA nanoball sequencing shows that the transcriptome is highly complex owing to numerous discontinuous transcription events. In addition to the canonical genomic and 9 subgenomic RNAs, SARS-CoV-2 produces transcripts encoding unknown ORFs with fusion, deletion, and/or frameshift. Using nanopore direct RNA sequencing, we further find at least 41 RNA modification sites on viral transcripts, with the most frequent motif, AAGAA. Modified RNAs have shorter poly(A) tails than unmodified RNAs, suggesting a link between the modification and the 3' tail. Functional investigation of the unknown transcripts and RNA modifications discovered in this study will open new directions to our understanding of the life cycle and pathogenicity of SARS-CoV-2.